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ABSTRACT 

Water is key in the evolution of protoplanetary disks and the formation of comets and icy /water 
planets. While high excitation water lines originating in the hot inner disk have been detected in 
several T Tauri stars (TTSs) , water vapor from the outer disk, where most of water ice reservoir is 
stored, was only reported in the closeby TTS TW Hya. We present spectrally resolved Herschel/YQFl 
observations of the young TTS DG Tau in the ortho- and para- water ground-state transitions at 557, 
1113 GHz. The lines show a narrow double-peaked profile, consistent with an origin in the outer disk, 
and are ~ 19 — 26 times brighter than in TW Hya. In contrast, CO and [C II] lines are dominated by 
emission from the envelope/outfiow, which makes H2O lines a unique tracer of the disk of DG Tau. 
Disk modeling with the thermo-chemical code ProDiMo indicates that the strong UV field, due to the 
young age and strong accretion of DG Tau, irradiates a disk upper layer at 10-90 AU from the star, 
heating it up to temperatures of 600 K and producing the observed bright water lines. The models 
suggest a disk mass of 0.015-0.1 M0, consistent with the estimated minimum mass of the solar nebula 
before planet formation, and a water reservoir of ^ 10^ — 10'^ Earth oceans in vapour, and ^ 100 times 
larger in the form of ice. Hence, this detection supports the scenario of ocean delivery on terrestrial 
planets by impact of icy bodies forming in the outer disk. 

Subject headings: astrochemistry - ISM: molecules - protoplanetary disks - stars: individual (DG Tau) 



1. INTRODUCTION 

Protoplanetary disks are the birthplaces of planets, 
thus the study of their physical and chemical structure 
is fundamental to comprehend the formation of our own 
solar system as well as of extra-solar planetary systems. 
One of the most intriguing issues on planet formation 
concerns the origin of oceans on Earth. It was argued 
that Earth formed "dry" and that ocean water was de- 
livered by impacts of icy bodies/protocomets originating 
from the cold outer disk, where most of the m ass (and 
water reservoir) is located (iMatsui fc Abel[T986[ ). To ad- 
dress this issue, many efforts have been devoted to ob- 
serve water in protoplanetary disks and to characterize 
its abundance and spatial distribution. 
In the hot dense inner disk region inside the so- 
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called 'snow line' where Tdust '~ 150 K, i.e. for radii 
smaller than ~ 1 — 3 AU in disks around T Tauri 
Stars (TTSs) (ILecar et al.l 120061) . ice cannot exist on 
dust grains and gas-phase chemistry converts all oxy- 
gen into water on timescales short compared to the disk 
evolution timescale. Beyond the snow-line, instead, wa- 
ter molecules will be frozen onto dust grains. However, 
(inter)stellar UV and X-ray radiation can penetrate in 
the disk upper layers and photo-d esorb a fraction of wa- 
ter ice back into th e gas phase (jCeccarelli et al.l [20051 : 
iDominik et al.|[2005[ ). The released water vapour may be 
eventually dissociated and re-formed in the gas-phase. 

H2O lines with upper level energies i?up > 1000 K, 
tracing hot water vapor in the inner disk regions, 
have now been observed in a number of protoplane- 
tary disks thanks to ground-based and Spitzer near- and 
mid-infrared obs e rvations (e.g. [Ca rr & Naiita 2008|; 
iSalvk et al.l 120081: iPontoppidan et aiTl2010a .bV. and, re- 
cently, far -infrared observations of the 63. 32 /im line with 
Herschel (jRiviere-Marichalar et al.ll2012t l. In contrast, 
cold water vapor at T < 200 K from the outer disk 
surface has turned out surprisingly difficult to detect 
in TTSs. Herschel/PACS detected the low excitation 
H2O 179.5 ^m hne {E^p = 114 K) only in jet-driving 
stars but due to the lack of spatial and velocity informa- 
tion, it is uncle ar if it originates in the disk or in the en- 
velope/outflow (jPodio et al.ll2012l ). Up to now, firm evi- 
dence for a cold disk water reservoir has been found only 
in the nearby (d - 50 pc) TTS TW Hya, through the 
detection of the fundamental ortho an d para lines at 557 
and 1113 GHz with Herschel/lllFl (jHogerheiide et al.l 
I20TI . While 0-H2O 557 GHz line profiles in Class 
and I sources show velocities of ~ 11 — 138 kms~^, and 
'--^ 5 — 54 kms~^, suggesting they are dominated by emis - 
sion from the envelope/outfiow (jKristensen et al.1l2012| ). 
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the H2O emission from TW Hya sliows a narrow single- 
peaked profile {FWHM - 0.96 - 1.2 kms~^) consistent 
with an origin in the face-on disk. A hidden reservoir of 
icy bodies of 1.5 M^, equivalent to several thousands of 
Earth oceansH is inferred. Additional studies are neces- 
sary to investigate this hypothesis, but only upper limits 
were obtained toward a couple of other TT Ss targeted 
with HIFI, e.g. DM Tau ()Bergin et al.ll2010[) . 

DG Tau is a young TTS at 140 pc associ- 
ated with particularly strong accretion/outflow activ- 
ity (e.g. .Hartigan ct al. 1995; Dougados ct al. 2000), 
and where we previously detected unresolved emission 
in the H2O 78.7, 179.5 /^m lines with Herschel/PACS 
(jPodio et al.l 120121) . In this Letter, we present clear de- 
tections of the H2O 557, 1113 GHz lines towards this 
source. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed DG Tau (aj2ooo = 04^ 27'° 04!7, (5j2ooo = 
-1-26° 06' 16'.'3) with the Heterodyne Instrument for the 
Far Infrared (HIFI, dc Graauw ct al. 2010) on board the 
Herschel Space ObservatorjQ (Pilb ratfet al.ll2010t ). The 
observations target the two fundamental water lines, o- 
H2O lio-loi and P-H2O lii-Ooo, the ^^CO (hereafter CO) 
and "CO 10-9, and the [C II] 2P3/2-^Pi/2 lines (OB- 
SID: 1342239630, 1342250208, 1342249594, 1342249646). 
They were acquired in the HIFI bands 1, 4, 5, and 
7, with a single on-source pointing and in dual beam 
switch mode with fast chopping 3' either side of the tar- 
get. The Wide Band Spectrometer (WBS) and the High 
Resolution Spectrometer (HRS) were used in parallel, 
with a spectral resolution of 1.10 and 0.25 MHz, respec- 
tively. The Half Power Beam Width (HPBW) ranges 
from 11" to ~ 38", depending on frequency. 

HIFI data were reduced using HIPE $3. Fits files 
from level 2 were then created and transformed into 
GILDA^3 format for data analysis. The spectra were 
baseline subtracted and then resampled at 0.6 kms~^ to 
increase the sensitivity. Note that the V-spectrum of the 
0-H2O and P-H2O lines are affected by ripples, degrading 
the quality of the baseline and resulting in an rms larger 
than the one measured in the H-spectrum. Therefore, 
in the following, we will analyse the 0-H2O and P-H2O 
emission based solely on the H-spectrum. 

The HIFI dataset is complemented by observations 
of the CO 3-2 line performed on January 2010 at the 
JCMT 15m telescope (Manua Kea, Hawaii, USA) using 
the HARP-B heterodyne array and ASCIS correlator, 
providing a spectral resolution of 0.25 kms^^. The spec- 
trum was resampled at 0.6 kms~^ to be compared with 
the HIFI data. 

Antenna temperatures, Ta, are converted to mean 
beam temperature, Tmh (for HIFI mean beam efficiency 
are bv lRoelfsema et al1l2012j ). Integrated line intensities, 

/ T^bdV, and line fluxes, F^bs = x / T„,bdV x 

11 1 M0= 5.971027 g, 1 Earth ocean ~ 1.5 lO^^ g 
^2 Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 

HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Sci- 
ence Center, and the HIFI, PACS and SPIRE consortia. 
" http://www.iram.fr/IRAMFR/GILDAS 



TABLE 1 
Lines integrated intensities 



Transition" 






HPBW 






GHz 




n 


K kms ^ 


0-H2O lio - loi 


556.936 


0.76 


38 


0.10 ± 0.01 


P-H2O 111 - Ooo 


1113.343 


0.74 


19 


0.12 ± 0.01 


CO 10-9 


1151.985 


0.64 


18 


5.8 ± 0.1 


"CO 10-9 


1101.350 


0.74 


19 


0.28 ± 0.01 


[C II] 2P3/2-2Pi/2 


1900.537 


0.69 


11 


3.1 ± 0.2 


CO 3-2 


345.795 


0.66 


14 


32.5 ± 0.4 



"All lines are observed with //ersc/ie//HIFI except CO 3—2 which 
is observed with JCMT/HARP-B 

^Frequenc ies are from the J et Propulsion Laboratory molecular 
database llPickett et al.lll998l) 




, are summarized in tables [T] and [3l 



3. RESULTS FROM OBSERVATIONS 

The observed line profiles are shown in figure [TJ The 
JCMT CO 3-2 line profile in panel (a) suggests that the 
systemic velocity is Kvs ^ -1-6.2 km s~\ consistent with 
previous studies (Schu ster et al.l Il993t iKitamura et al.l 
[1996; Testi ct al. 2002). 

We detect both the ortho-H20 lio-loi 557 GHz and 
the para-H20 In-Ooo 1113 GHz lines (Kp ^ 61, 53 K) 
with a signal-to- noise of 10 and 12, respectively. They 
are centered at the systemic velocity and show a narrow 
double-peaked profile (FWHM ~ 5 — 6 kms^^). 

CO 10-9, "CO 10-9, and [C II] 158 ^lm lines have a dif- 
ferent profile than H2O, with a single peak near systemic 
velocity (at t^LSR — +6.5 kms^^ in CO and -1-5.5 kms~-^ 
in [C II]), and a pronounced blue wing extending down 
to and —5 kms^^ respectively (i.e. 6 and 11 kms~^ 
away from systemic). The bulk of CO and [C II] emis- 
sion close to systemic velocity likely originates in the en- 
velope, as suggest ed by the ^'^CO 2-1 channel maps by 
iTesti et all (|2002f ) that indicate emission extended over 
~ 10" at velocities — Kiys|< 1-5 kms~^. The observed 
blue wing, instead, may originate in a slow outfiow, per- 
haps linked to envelop e dispersal motions as proposed by 
IKitamura et"all (fT996h . For the [C II] 158 /im line an ori- 
gin in an extended structure is further confirmed by the 
fact that the fiux in the HIFI beam of ^ 11" is ^ 4 times 
lower than the total co-add ed flux in the Herschel/PACS 
observations (47"x47". iPodio et al.|[20ia V 

On the other hand, several arguments suggest that 
the H2O emission is compact and is likely dominated 
by emission from the outer region of the protoplanetary 
disk of DG Tau, and not from the envelope/outflow: 

• The H2O line profiles are different from those of 
CO and [C II] observed with single-dish telescopes. 
They are much more symmetric about the systemic 
velocity and do not show the extended blue wing 
seen in these other tracers. 

• The peaks of the H2O line profiles coincide with 
the two narrow velocity ranges (J^ — Kiysl= 1-5 — 
2.5 kms^^) where ^"^CO 2-1 interferometric maps 
show compact emission with a velocity gradient 
perpendic ular to the jet ax is, consistent with disk 
rotation (|Testi et al.ri2002f) . The "CO 2-1 line 
profile obtain ed by integrating the interferometric 
maps bv .Testi et all (|200a) over a 2" beam, i.e. by 
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Fig. 1. — HIFI spectra of (a) 0-H2O lio-loi (grey histogram); 
(b) P-H2O lii-Ooo (grey histogram); (c) CO 10-9, "CO 10-9, and 
[C II] ■^P3/2-'^Pi/2 (green, red, and blue histograms, respectively). 
In panel (a) JCMT CO 3-2 (red histogram), and "CO 2-1 pro file 
obtained integrating interferometric maps bv l'l'esti et all 120021) on 
a 2" beam (blue histogram) are also shown. The vertical dotted 
line indicates the systemic velocity (Vlsr = -1-6.2 kms~^). The 
H2O line profiles predicted by the "low dust opacity" ProDiMo 
disk model are overplotted (magenta lines). The 0-H2O line flux 
is underpredicted by the model by a factor ~2.2, hence the line 
proflle is multiplied by this factor to help the comparison with 
observations. 

cutting out any extended component, is similar to 
the H2O line profiles, with peaks at the same ve- 
locities. In contrast, the CO 3-2 profile, obtained 
with the JCMT collecting all the emission in the 
14" beam, does not peak at the same velocity as 
the H2O and ^^CO compact component. This is 
particularly clear in the blue part of the profile. 

• Assuming keplerian rotation, and an inclination of 
i ~ 38° from the line of sight (Eisloffcl &: Mundt 



I1998D the peak separation of the H2O lines 
(AT^cp ^ 3 — 3.5 kms~^) indicates an outer disk ra- 
dius i?out(H20)- 77 - 105(M^/0.7Mp) AU. For a 
stellar mass of ^ 0.7 M0, as assumed in lTesti et al.l 
(l200l . the inferred i?out(H2 0) is in agreement 
with the disk outer radius, ^ 72 — 89 AU, esti- 
mated from sub-arcsecond dust continuum ni aps at 
1.3, 2.8 mm with CARMA (|Isella et al.ll20Tol ). The 
maximum velocities covered by the line profiles, in- 
stead, set an upper limit to the inner radius of the 
line emitting region i?in(H2 0)< 19 AU, since more 
extended line wings could be hidden in the noise. 

• The H2O line profiles are reproduced by an op- 
tically thick, vertically isothermal keplerian disk 
with Tox oc r^'^-^ viewed at 38° with an excitation 
temperature at i?out of 70 an d 32 K for the ortho 
and para lines respecti vely (jBeckwith &: Sargentl 
[1993I: rCabrit et al.l[2006l) . 

Given the evidence listed above, the fundamental water 
lines, even when observed with a 38"-19" beam, appear 
to be dominated by compact emission. Although we can- 
not exclude contamination from the outflow, which could 
explain the larger FWHM and the asymmetry of the o- 
H2O 557 GHz, the detected double-peaked H2O lines 
prove to be a good tracer of the outer protoplanetary disk 
of DG Tau, with less confusion from envelope/outflow 
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DG Tau shows emission a lso in high-excitat ion H2O 
hues observed with PACS (|Podio et al.l 120 12| ). With 
Eup ~ 200 — 1070 K these are thought to originate in an 
intermediate disk region between a few and a few tens o f 
AU from the star (e.g. iRiviere-Marichalar et al.ll2012D . 
The exception is the low-excitation H2O 179.5 fim line 
{Eup 114 K) which according to previous disk model- 
ing is predicted to form in the outer disk like the 557, 
1113 GHz lines (Kamp et al. submitted). The observed 
H2O 179.5 fim I 557 GHz line ratio is i?lobs = 25 ± 6, 
consistent with LTE optically thick emission in the 
Rayleigh- Jeans limit, i.e. for temperatures larger than 
a few hundreds K (i?LTE-thick ~ 27). On the other 
hand, the line ratio between the para- and the ortho- 
fundamental lines (i?2obs = 2.5 ± 0.3) is around three 
times lower than i?LTE-thick ^ 8. This can be explained 
if the lines are excited in a region where the gas density 
is lower than the lines critical density (~ 2 10^ and 
- 210** at 50 K for the 557 and 1113 GHz lines) 
and/or where the temperature is below their upper level 
energies. Also, the observed line ratio could be affected 
by emission from the envelope/outflow. 



4. MODELING H2O IN THE DISK OF DG TAU 

Detailed disk modeling is required to test the disk hy- 
pothesis and to derive an estimate of the water mass. 
The latter cannot be inferred from observations since the 
lines are likely optically thick. We include in our anal- 
ysis the fluxes and upper limits obtained for the water 
lines falling between 63.3 and 180.5 /im observed with 
PACS as part of the Hersch d Key Project GASPS (PI: 
B. Dent) (jPodio et al.ll2012D . The two detected 0-H2O 
lines at 78.7 and 179.5 /xm are spectrally and spatiall 
unresolved, thus their origin is unclear. In iPodio et al 
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Fig. 2. — Disk region from which 50% of the 0-H2O 179.5 fim (in 
green), 0-H2O 538.3 (or 557 GHz, in red), and P-H2O 269 fim 
(or 1113 GHz, in orange) Hne emission arises according to the "low 
dust opacity" disk model. The grey colour indicates the water 
density, njjjO (cm~^), the dotted black and blue curves the gas 
temperature and density, and the white solid curve the snow-line 
(i.e. Td^^t = 150 K). 

(HoTl) a shock -origin was favoured based on the large hne 
fluxes, which are difflcuh to reproduce with disk models 
for typical TTSs parameters. However, since the pro- 
files of the ground-state water lines are consistent with a 
disk origin, we test the predictions of a dedicated model 
for DG Tau by comparing them with observed H2O line 
fluxes and profiles. 

We use a parametrized disk model calculated with 
the thermo-chemical disk modeling code ProDiMo 
(jWoitke et all 120091: iKamp et all 12010( 1 . We adopt the 
stellar spectral type K7 {Tcs — 4 000 K) and veiling - 
corrected stellar radius of 1.8 Rq (|Fischer et al.l [20111 ). 
The resulting stellar luminosity ~ 1 Lq yields a stellar 
mass Mi, ~ 0.7 and an age of 2.5 1 0^ yr using the evo- 
lutionary tracks of lSiess et all (I2000D. To reprod uce the 
lUE UV/optical spectrum (jGuUbring et"alll2000f ). we set 
the UV excess fraction /uv = ^(910 - 2500 A)/L* = 0.2 
and adopt a power law slope L\ ss A"'' '^. We also 
account for the effect of X-ray r adiation from the stel- 
lar cor ona (Lx = 10'^° e r g s-MGiidel et "all 120071) fol- 
lowing [Areiuetin (|2011[ ): iMeiierink et all ()2012[ ). The 
disk inner and oute r radius are set to Rm = 0.16 AU 
()Akeson et alJ[200l . and i?out = 100 AU, in agreement 
with i?out(H20) inferred from the observed H2O pro- 
files. We assume the dust size distribution and disk 
dust mass from the "low dust opaci ty model" — a 50/50 
mixture of astronomical silicates (iDraine fc Led [1981 
and amorphous ca rbon (jZubko et al.l 119961 ) — used by 
llsella et all (|2010l ) to reproduce the observed 1.3, 2.8 
mm emission {n{a) « a""^ with q — 3.5, where a is 
the dust grain radius, the minimum/maximum grain size 
are amin = 0.005 /im and Omax = 5 cm). Using the 
standard dust-to-gas ratio of 0.01 the gas mass is set 
to 0.1 Mq. The disk is thought to be perpendicular 



to the jet, thus i = 38° (|Eisloffel fc Mundtl [19981) . We 
assume a parametrized disk shape with surface density 
E w r-\ and scale height = 0.008 AU (r/0.16 AU)^ ^ 
No dust settling is invoked, i.e. dust and gas are well 
mixed throughout this young disk. The PAH fraction 
is 0.01 with respect to the ISM abundance of 10^^-^^ 
PAH particles/H-nucleus. All parameters adopted for 
the model are summarized in table [2] 

The line profiles and fluxes are obtained by first solving 
the statistical equilibrium with 2D escape probability to 
obtain the level populations, and then using 2D radia- 
tive transfer (collision rates as listed in table 3 of Kamp 
et al., submitted). The region from which 50% of the 
H2O line emission arises, instead, is obtained using ver- 
tical escape probability and without accounting for disk 
inclination. Figure [2] indicates that the H2O 179.5 /^m 
line observed with PACS originates in the same region as 
the fundamental water lines at 557, 1113 GHz observed 
with HIFI, i.e. in an upper disk layer {z/r ^ 0.35 — 0.6) 
located at '--^ 10 — 90 AU distance from the star. In this 
region the gas temperature is ~ 50 — 600 K and water 
is formed mainly through gas-phase reactions and par- 
tially dissociated by UV photons and collisions with C+ 
and H+ . Including self- shielding for all photodissociating 
species produces at most 20% lower fluxes. The gas den- 
sity is 10^ — 10^" cm~'^ thus, as suggested by the observed 
H2O 179.5 /xm/557 GHz line ratio, these lines are close to 
LTE and optically thick (r ~ 10^ - lO"*). The ortho-to- 
para ratio (OPR) is calculated from the gas temperature 
at thermal equilibrium and is 1.5-3 in the line emitting 
region. However, since the H2O lines are optically thick, 
the model results are not dependent on the OPR. 

As shown in figure [1] the model reproduces the 
P-H2O line flux and proflle, and the ratio o- 
H2O 179.5 /im / 557 GHz is i?lmod — 25, in agreement 
with the observed value. On the other hand, the ob- 
served ortho lines at 179.5 /im and 557 GHz are underpre- 
dicted by a factor 2. As a consequence the observed p- 
H2O 1113 / 0-H2O 557 line ratio is overpredicted by a fac- 
tor 2.4 (i?2,„od = 6.1). Kamp et al. (submitted) discuss 
in detail the uncertainties when modeling water emis- 
sion in disks. They show that the assumed surface chem- 
istry, adsorption energy and photo-desorption yields, and 
metal abundances can affect water line fluxes by a factor 
of a few. In particular, the low-excitation water lines are 
very sensitive to the adopted radiative transfer method 
and to the uncertainties in the collision rates. Moreover 
the disk m odel is not a.ccoun ting for the X-ray emission 
by the jet (jGiidel et alJl2008t ) which illuminates the disk 
surface from above. This may boost water formation 
through H^O^ recom bination {H^O^ +e~ ^ H2O + 
IMeiierink et alj[2012i) . In general, the model can repro- 
duce all the H2O lines observed with PACS and HIFI 
within a factor 2 (see table [3]). The emission in the CO 
and [C II] lines is predicted to be a factor 3-9 lower than 
observed, which is in agreement with the observed pro- 
files indicating that the bulk of the emission originate 
from the envelope/outflow. 

The disk model indicates that the disk contains ~ 0.4 
Mg of water vapour, and two orders of magnitude larger 
mass in ice: M(H2 0#) ~ 100 M®. To understand the 
reliability of the estimated water mass in the disk, we cal- 
culate a second model assuming the dust size distribution 
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TABLE 2 



IjUW U U b i (Jr^i-Vt^i i 1 UioiV MUUiljij. 


STAR AND DISK PARAMETERS 


Effective temperature 


Tc« (K) 


4000 


Stellar mass 


M, (Mq) 


0.7 


Stellar luminosity 


L, (Lq) 


1 


UV excess 


/uv 


0.2 


UV power law index 


Puv 


-0.3 


X-rays luminosity 


Lx (ergs~l) 


1 n30 
10 


Disk inner radius 


Rin (AU) 


0.16 


Disk outer radius 


Rout (AU) 


100 


Disk dust mass 


Mdust (Mq) 


1 10" 


Dust-to-gas ratio 


dust-to-gas 


0.01 


Solid material mass density 


Pdust (g cm-3) 


3.5 


Minimum grain size 


Imin (^m) 


0.005 


Maximum grain size 


dmax (cm) 


5 


Dust size distribution index 


1 


3.5 


Disk inclination 


i n 


38 


Surface density E r~'^ 




-1 


Scale height at 


Ho (AU) 


0.008 


Disk flaring index H{r) = Hq 


1.2 


Fraction of PAHs w.r.t. ISM 


/pah 


0.01 



The total water reservoir, M(H2 0)gas-i-ice ^ 7 — 100 M^, 
is a factor of a few up to two orders of mag nitude larger 
than for TW Hya (jHogerheiide et al.|[20lll) . 



TABLE 3 

Observed and disk-model predicted H2O fluxes 



Line 


A 






-^mod 




fim 


K 


Wm-2 


Wm- 


2 


PACS Observations 


0-H2O 


63.3 


1070 


< 4 10-^' 


1.9 10" 


17 


0-H2O 


71.9 


843 


< 1 10-1'^ 


1.7 10- 


17 


0-H2O 


78.7 


432 


1.9 ± 1.4 10-1'^ 


2.1 10- 


17 


0-H2O 


179.5 


114 


1.5 ± 0.3 10-" 


7.7 10- 


18 


0-H2O 


180.5 


194 


< 1 10"" 


3.2 10- 


18 


P-H2O 


78.9 


781 


< 1 10"" 


9.7 10- 


18 


P-H2O 


89.9 


297 


< 1 10-1'^ 


1.4 10- 


17 


P-H2O 


144.5 


396 


< 1 10"" 


2.7 10" 


18 
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an d disk dust mass fr om the "high dust opacity model" 
by llsella et alj (j2010D . This implies around an order of 
magnitude lower dust mass in the disk and consequently 
around an order of magnitude lower gas mass, and water 
vapor and ice mass (Mgas = 0.015 M©, M(H20)-- 0.06 
Mq, M(H20#)-- 7 Me). We find that this model can 
reproduce equally well the observed H2O line fiuxes, be- 
cause the "high dust opacity model" implies lower opac- 
ity at UV wavelengths and thus a deeper UV penetration 
in the outer disk regions. Hence, the dust size distribu- 
tion is crucial to constrain the disk mass and water reser- 
voir, leading to an uncertainty of one order of magnitude. 



5. CONCLUSIONS 

The present detection of the 0-H2O and P-H2O lines at 
557, 1113 GHz in the TTS DG Tau is crucial for several 
reasons: (i) so far, emission in the fundamental water 
lines has been observed only in one TTS, TW Hya; (ii) 
we detect for the first time a double-peaked profile in the 
H2O lines, which is a strong kinematic evidence for an 
origin in the outer disk (from ^ 10 — 90 AU); (iii) wa- 
ter turns out to be a unique tracer of the protoplanetary 
disk of DG Tau, because is less contaminated by enve- 
lope/outflow emission than CO lines; (iv) once corrected 
for distance the H2O lines are ~ 19 — 26 times brighter 
than in TW Hya. According to our models, the reason is 
the 10 times higher UV flux of DG Tau, which heats the 
outer disk surface layer up to temperatures of ^ 600 K 
(only ^ 30 K in the case of TW Hya). In addition, 
the disk around DG Tau is more massive and compact 
leading to higher volume densities in the surface layers, 
which makes the warm neutral chemistry even more ef- 
ficient; (v) the adopted models suggest a disk mass of 

0. 015-0.1 M0, depending on the assumed dust size dis- 
tribution, and a water reservoir (gas-|-ice) of 7-100 M®, 

1. e. at least a factor of a fe w larg er than estimated for 
TW Hya (|Hogerheiide et aLllMTl ). 

While the inferred disk mass is consistent with the min- 
imum mass of the solar nebula to form our solar system, 
the detection of water vapor in the outer region of the 
disk, where comets are believed to form, and the esti- 
mated water mass of a few ~ 10* — 10^ Earth oceans, 
supports the scenario of impact delivery of water on ter- 
restrial planets by means of icy bodies. 
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